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Abstract 

The resummation of threshold logarithms for direct photon pro- 
duction cross sections in hadronic collisions is presented. The resum- 
mation is based on the factorization properties of the cross section and 
is formulated at next-to-leading logarithmic or higher accuracy. Full 
analytical and numerical results for the next-to-next-to-leading order 
expansion of the resummed cross section are given. A substantial 
reduction of factorization scale dependence is observed. A compari- 
son to experimental results from the E-706 and UA-6 experiments is 
presented. 
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1 Introduction 

The production of photons with large transverse momenta has long been 
recognized as an important probe of hard-scattering dynamics. In particular, 
the subprocess qg — > jq offers sensitivity to the gluon distribution function 
in lowest order. Thus, direct photon production is often cited as having the 
potential to place strong constraints on the gluon distribution in global fits 
of parton distributions. In recent years, however, the situation has become 
more complex. Recent detailed phenomenological studies have shown that 
it is not possible to describe the existing set of direct photon experimental 
results using a next-to-leading order QCD calculation JTJ, [J. Even with the 
full freedom of separately adjusting the renormalization and the factorization 
scales for the distribution and fragmentation functions, one can not correctly 
describe all of the experimental results . In 0] it was pointed out that for 
many of the data sets the theoretical predictions fall below the data near the 
low-pr end of the pt range spanned by the data. This has been interpreted [|J 
as providing evidence for the inclusion of fc^-smearing which is due to initial 
state radiation. However, this interpretation is not completely supported 
by all of the available data and the analysis of Jl] shows that agreement 
with a subset of the data can be obtained if one considers a restricted set of 
data taken only on proton targets (so as to remove possible nuclear effects) 
and also if minimum px cuts on the data sets are applied which restrict the 
comparisons to regions where the scale dependence of the theory is moderate, 
and, hence, the theoretical predictions are thought to be reliable. 

In light of the situation described in the preceeding paragraph, efforts are 
underway to improve the precision of the theoretical predictions for direct 
photon production. One such type of improvement involves the resummation 
of large logarithms which result from phase space limitations on real gluon 
emission near the threshold region for the underlying parton scattering sub- 
processes. Such corrections are expected to be important in the region of 
large photon transverse momentum, where xt = 2px/VS approaches one. 
Nevertheless, it is of interest to examine the effects of such resummation 
calculations over an extended range of Xt- 

Recently, threshold resummation for direct photon production has been 
discussed in Refs. || |6], [7|, §. In Refs. [|, [Tj, a resummation calculation was 
presented for the inclusive Et distribution, da/dEx, and compared with some 
of the available fixed target data. The results showed that, as expected, the 
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corrections are large as xt approaches the edge of phase space. In addition, 
the scale dependence of the resulting distributions is reduced compared to 
that of the next-to-leading order calculation. Finally, in the region containing 
the data the resummed cross section was nearly the same as the next-to- 
leading order cross section if the factorization and renormalization scales 
were chosen to be \i = Pt/2- 

Another formalism exists for treating the resummation of threshold loga- 
rithms in inclusive processes, a review of which can be found in Ref. ||. This 
formalism has been applied to lepton pair, heavy quark, and dijet production 
at fixed invariant mass and more recently has been extended to single-particle 
inclusive cross sections [§], including direct photon || and W + jet produc- 
tion ITTjL ITTJ. In the present paper this formalism is applied to the inclusive 



invariant cross section for direct photon production. One important feature 
of this method is that one directly obtains the angular dependence (or ra- 
pidity dependence) for the produced photon. The resummation explicitly 
includes all the leading and next-to- leading logarithms (NLL). For detailed 
numerical work and comparison to other calculations, the resummed cross 
section expression is expanded to next-to-leading order (NLO) and next-to- 
next-to-leading order (NNLO). Matching the resulting expressions to existing 
complete NLO calculations yields the leading, next-to-leading, and next-to- 
next-to-leading logarithms (NNLL) through 0(aa^). The NLO expansion 
agrees fully near partonic threshold with existing exact NLO calculations 
1^, [13| while our NNLO results provide new predictions. A comparison with 



recent high statistics fixed target results is also presented. 



2 Factorization and resummed cross section 
2.1 Factorized cross section 

In this section we discuss the application of the threshold resummation for- 
malism described in Ref. || to direct photon production. The self-contained 
presentation follows the techniques of Refs. [|5], |TU|] . 

The quantity to be calculated is the inclusive invariant cross section for 
the production of photons in hadronic collisions: 



hAipA) + h B {p B ) -> 7(p 7 ) + X 



(2.1) 
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Fig. 1. (a) Factorization for direct photon production near partonic thresh- 
old, (b) The soft-gluon function S, in which the vertices c link ordered 
exponentials, representing the partons in the hard scattering. 



The Mandelstam invariants formed form the hadron and photon four- vectors 
are 

S=(pa+Pb) 2 , T=(p A -p 7 ) 2 , U=(p B -p 1 ) 2 . (2.2) 

The parton-parton scattering subprocesses contributing to direct photon pro- 
duction in lowest order are 



q(pa) + g(Pb) — ► 7(p 7 ) + q(pj) 



(2.3) 



and 



q(Pa) + q(Pb) — > 7(P 7 ) + 9(pj) ■ (2-4) 
The corresponding Mandelstam invariants constructed from the parton and 
photon four-vectors are given by 



S=(Pa+Pb) 2 , t=(p a ~p~ / ) 2 , U={p h ~p 1 f 



(2.5) 



which satisfy s i = s + t + u = 0dX threshold. The variable s 4 is the square 
of the invariant mass of the system recoiling against the observed photon. It 
will also be useful to define a corresponding variable, S4 — S + T + U, using 
the hadronic Mandelstam variables. 

The factorized form of the cross section for direct photon production is 
a convolution of the parton distribution functions <f) with the parton-parton 
hard scattering cross section: 



E, 



= J dx a dx b <p fa/hA (x a , /4) <p h/hB (x b , n%) 



5 



x E f a ^ lX {s, t, u, fx F , a fl (/4)) . (2.6) 

The initial-state collinear singularities have been factorized and absorbed 
into the 0's at a factorization scale /ip, while is the renormalization scale. 

The threshold for the parton-parton scattering subprocesses occurs at 
S4 = 0. At this point the photon recoils against a massless parton and 
receives the maximum amount (one half) of the parton center of mass energy 
allowed. The values of x a and Xb corresponding to this point are the minimum 
values which can give rise to a photon of the specified rapidity and transverse 
momentum. Near threshold the phase space for the emission of real gluons is 
limited and there is an incomplete cancellation of infrared divergences against 
virtual gluon emission contributions. Although the remainder is finite, there 
are large logarithmic corrections. In general, a includes plus distributions 
with respect to s 4 at nth order in a s of the type 



ln m (s 4 A4) 
s 4 



m < 2n - 1 , (2.7) 



where Pt = (tu/s) 1 ^ 2 is the transverse momentum of the photon. 

To organize these plus distributions we introduce a refactorization in 
terms of new parton distributions ip and a jet function J that absorb the 
collinear singularities from the incoming partons and outgoing jet || [14], [15 



respectively; a soft gluon function S that describes noncollinear soft gluon 
emission JTJ], DJ; and a short-distance hard scattering function H. This 
refactorization is shown in Fig. 1 for the partonic cross section with the 
colliding hadrons in Eq. ( [2.6|) replaced by partons f a , fb- In our soft-gluon 
approximation the partons in the hard scattering are represented by eikonal 
lines (ordered exponentials). The cut diagram in Fig. 1 shows contributions 
from the amplitude and its complex conjugate, with H = hh*. 

The momenta of the parton lines entering the hard scattering subprocess 
are then x a pA and XbPB, where now pa and ps denote the momenta of external 
partons. If soft radiation of total momentum k$ is emitted, then momentum 
conservation gives x a pA + XbPs = P-y + Pj + k s where pj represents the 
momentum of everything in the final state excluding the photon and the 
soft radiation. Using the definitions of the Mandelstam variables presented 
previously, and neglecting corrections of order Sf , the relation between S4 
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and S4 can be expressed as [p] 

5 4 1 r 
~S ~ S 



%Pj ■ k s + p] + (1 - x a )2p A -pj+(l- x b )2p B ■ pj 



(-u) (-t) 

= W s + Wj + W a \-w b 

s s 

= fi-(l-x a )--(l-x 6 )-, (2.8) 
s s s 

where the dimensionless w^s measure the fractional contributions to 5*4/5 of 
the gluon emission associated with the functions ipi, J, and S. The refactor- 
ized partonic cross section can be written as a convolution || 

E~ d af j/p~* jX = H J dw a dw b dwjdw s i>f a /f a (w a )'ipf b/ f b (w b ) J(wj) 

x S{w S \/s/fi F ) 5 ( ^ -w 8 -wj- - w b ^^- j . (2.9) 

\ b s s J 

If we take moments of the above equation, with N the moment variable, we 
can write the partonic cross section as 

/ ^e- NS ^ E / af ^ X = ^ fa/fa (N a )^ fb/fb (N b )J(N)HS(V-s/m F )), 

(2.10) 

where 

$f a /f«(N a ) = J dw a e- N ^ fa/fa (w a ) , (2.11) 

and similarly for the other functions, with N a = N(—u/s) and N b = N(—t/s) 
for either partonic subprocess. 

Next, we take moments of Eq. ( |2.6| ) with the incoming hadrons replaced 
by partons and, using the last line of Eq. (|2.8|), we find 



= h/U mh/A (N b )E, d3&faf ^ X (N) ■ (2-12) 
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Comparing Eqs. (|2.12|) and ( |2.10| ) we may then solve for the moments of the 
perturbative cross section a: 



E. 



d^f afb ^x{N) _ j> fa , U {N a )j> fb/fb {N b ) ~ 



d 3 p^ 



ha/fa(Na) h b /f b (N b ) 



J(N)HS(y/s/(N(j, F )). (2.13) 



The moments of the plus distributions with respect to S4 in a give powers of 
IniV as high as ln 2n iV. In the next section these logarithms are resummed 
to all orders in perturbation theory. 



2.2 Resummed cross section 

The resummation of the iV-dependence of each of the functions in the refac- 



torized cross section, Eq. ( 2.13 ), is based on their renormalization properties 
and has been extensively discussed in Refs. |5], [LT], O, |T5|, The resummed 
cross section in the MS factorization scheme is given by: 



E, 



d 3 p~ 



exp £ [E^\N t \ 



2 r*Pi<d£ 



exp{^(j)(iV)} H (s,t,u,a s (/i F 



7/ l (« s (^ 2 ))-7/ 1 / I (^,«s(/i' 2 )) 



~ / \ rVs/N du! / 

x S (a s (s/N 2 )) exp ^ 2 R e r 5 (a e (ji 



(2.14) 



where £ M = Pj/y/s. The first exponent iv^(iVj), which resums the N- 
dependence of the ratio ip/4>, is given by |T^, [L7|, [L8| 

r N t -i _ l 



E^(Ni) 



dz 

1 — 2 



+ -v 



C/0 



(1-z) 2 A 
a s {{l-zf{2 Pl .Cf) 



AW [a, (\(2 Pl ■ C) 2 ) 



where, at next-to-leading order accuracy in In N, we have 



a.. 



7T 



7T 



(2.15) 



(2.16) 
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and 



7T 



(2.17) 



Here Cf = Cp= (JV| — l)/(2iV c ) for an incoming quark, and C/ = Ca = iV c 
for an incoming gluon, with N c the number of colors, while 



K = C A 



'67 _tt^ 
18 6 
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(2.18) 



where n/ is the number of quark flavors [[19 



The anomalous dimensions 7/ and 7// of the fields ip and (T^, 0, ^0 
are given by 



7? 



3 a s 
-Gf — 

4 7T 
4 7T ' 



7w = - In i\T-- C>^, 



a. 



7T 



7gg = — ( Ca hiN — 



7T 



(2.19) 



for quark and gluon jets, respectively, where (3q = (11CU — 2n/)/3 is the 
one-loop coefficient of the j3 function. 

The exponent Ejn, which resums the iV-dependence of the final-state jet, 
is given by || [10], [15 



E[j)(N) 



1 z n-i _ 1 



dz- 



1-z 



(i _ z)2 y A (J) [a s (A S )] 



+ Bfjj [a.((l - *)*)] + B'{ J} [a.((l - z) 2 *)] } , (2.20) 
with given by Eq. ( |2.16| ) and -B(j) given for quarks by P, PHI 

3^ 



B[ q) = ^(-1)C F , B» q) = ^C F iM2v q )-l] 
and for gluons by 



a 



0o 



a. 



(y) 



7T 



7T 



The z/j terms are gauge dependent. They are defined by 

(A-n) 2 



Ui = 



n 



(2.21) 



(2.22) 



(2.23) 



where = pi^2/s are the particle velocities and n is the axial gauge vector, 
chosen so that Pi ■ ( = p, t ■ n for i = a, b ]5|, [TO . 
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Fig. 2. One-loop eikonal vertex corrections for partonic subprocesses in 
direct photon production; the eikonal lines represent the incoming partons 
and the outgoing jet. 



2.3 Soft anomalous dimensions 



The evolution of the soft function in Eq. 
renormalization group equation 



(p.!4|) , which follows from the 



' d d" 



S = -2(ReT s )S. 



(2.24) 



is given in terms of Ts, the soft anomalous dimension. is calculated 
explicitly at one loop by evaluating the eikonal vertex corrections in Fig. 2. 
The color basis for the hard scattering consists of only one tensor, c = Tp, 
and there are three eikonal lines connecting at the color vertex. Note that the 
color structure of the hard scattering for direct photon production is much 
simpler than for heavy quark or dijet production || |TJ], [TS], [H], |22|; hence, 
here Ts is simply a 1 x 1 matrix in color space. 

The soft anomalous dimension for the process q(p a )+g(Pb) l(p-y)+q(pj) 
is given by |5], ||, |T0|1 



1 s 



'79 



2tt 



2 In 



In 



— u 

s 
t 

u 



m(2z/ a 



1 



Til 



The soft anomalous dimension for the process q(p a )+q(Pb) 
is given by §, |, 

{C F [- ln{4v q Vq) + 2- 2m] 



l(Pi 



(2.25) 

^g(pj) 



1 S 



2tt 
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(2.26) 



Note that the soft anomalous dimensions are complex and that they also 
include the gauge dependent z/j terms; however, all gauge dependence cancels 
out explicitly in the resummed cross section. 

Finally, note that the resummed cross section, Eq. ( |2.14| ), can be rewrit- 



ten in a form |T(| more convenient for the calculation of the fixed-order 
expansions, as 



E, 



d 3 Vf a f b ^x(N) 
d 3 p^ 



x 



EW(Ni) - 

x exp{E[ J) (N)} exp 
x s(a s (s/N 2 )) exp 



H [a s (ji 2 R ) J exp 



2 i^^pm^)) 

Mi? f 1 



i z n-i _ j 
dz 

1 — 2 



2Rer 5 (a s ((l-z) 2 s) 



(2.27) 



3 NNLO expansion of the resummed cross 
section 

The formalism presented in the preceeding section results in expressions for 
the moments of the resummed cross section. An inverse Mellin transform is 
needed in order to obtain results for the pt and y dependence of the cross 
section. Such a transform involves integrals over regions of phase space where 
the argument of the running coupling vanishes and one must adopt a pre- 
scription for treating the resulting singularities; different prescriptions may 
give different results. On the other hand, it is possible to expand the mo- 
ment expressions in powers of as and perform the inversion order-by-order 
at fixed as- In this way one avoids having to adopt a specific prescription for 
the inversion. Of course, one does not have the fully resummed result after 
such a procedure. However, the work of Ref. shows that the expansion 
converges quickly until one reaches the region very near xt — 1. Therefore, 
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in regions relevant to current experiments, the corrections beyond NNLO 
are expected to be small. At any rate, this is a useful exercise for several 
reasons. First, one may check the resulting NLO 0(aa 2 s ) expressions with 
the corresponding terms from existing complete NLO calculations, thereby 
obtaining a test of the formalism. Second, one may generate analytical and 
numerical predictions for the dominant terms near threshold at NNLO. These 
terms represent the leading, next-to-leading, and next-to-next-to-leading log- 
arithms which would be present in a full 0{aa z s ) calculation. 
The results can be cast in the common form given by 



d 3 a MS 



ij^l+X B 



X 




= (Tij (S, t, U) 

2k-l 



4s(s^) + J2 c i+2 



1=0 




(3.1) 



where, for ease of notation, the subprocess labels have been suppressed on 
the c coefficients. It should be noted that these coefficients depend in general 
on s, t, and u. The expressions for the two Born terms are given by 



B 1 aa s 2 / s t 



and 

a * = ^;— e Au + v (3 - 3) 



where e q is the charge of a quark of type q. The 0(aa 2 s ) results for the 
coefficients c are given by 



cl(qg) = C F + 2C A 
c\{qg) = ~C F - (C F + C A ) In (j| 



c\(qg) = [-^-hj F -C F \n[-j)-C A \n 





h- §m(f ] (3.4) 
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and 



c\(qq) 

c\(qq) 



4C F - C A 



'/4 N 
pV 



C F -- -l n li In qi + q ln ^ 



.Pt. 



(3.5) 



These expressions are equivalent to the corresponding terms in the exact NLO 
results quoted in Ref. |TB[. The comparison is most easily made by rewriting 
the results in terms of the variables v and w used in |13| as is shown in the 
Appendix. Note that, with respect to the 5(s^) contribution, the expansion 
reproduces only the scale-dependent terms in c\. The full 5(s^) terms can be 
simply read off from the exact NLO calculation. 

The NNLO corrections in the MS scheme for the qg — > ^q subprocess are 
given by the following coefficients: 



4(q9) 

cl(qg) 



with 



and 



l - (C F + 2C A f 



--(C F + 2C A ) -C F + (C F + C A ) ln ^ 



4(qg) = (C F + 2C A ) 



x 4 

c T m (|) +c r 



.Pt, 



4 



+ 



~C F + (C F + C A ) In ^ 



2 \l2 



+ (C F + 2C A )^h, 



A, 



+ \k{C f + 2C A ) + ^-(3 C F 
2 lb 



799 



2C A ? 



Pa 3 

— -^-C^ln 



-C. ln 



-u 



.Pt , 

(3.6) 
(3.7) 



jqg 



/4 



I 11: 1 _ lev 



--CVln 2 



-N c In 2 

4 c 



+ ^C F ln 



s + t 



-C F ln 
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+ 



T, 



'19 



s 



t - 

+ - In 

s 



s + t 

s 

s + t 



t (2s + t 



T, 
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+ 2s 
1 (Zs + i 
2 



+ 



In 



7T 2 



-t 

s 



In 



2s + t 
s 

s + t 



s 
2t 



7T 



In — + — 1 - 4 



t (2s + t 
4s~ 



s + t) (s + tf 



In 2 



-t 

s 



+ 5 



+ 2 



(s + t) 



s + t 



t 



+ ^ In' — - 



-t 



+ l(^ 



(s + t) 2 t 2 



2t 



In 2 



+ - I In ( — ) In 



s + t 
s 

s + t 



(3.8) 



The corresponding coefficients for the qq — > 75 subprocess are: 



c\(qq) 

+ (AC F 
with 



- (ac f - c A y 



-3 (AC F - C A ) C F In - I (5CV - \C 



(AC F - C A ) 



c'f In 



+ 2C F In 



.A) 



/4 



7T 



l + («7,- C ,) T + S - T («7, 



/4 



C^) 2 (3.9) 



C, 



2 1 s 



(3.10) 



and 



-f 



2 2 2 



—u 



2 I Pt 



6 ) 



p¥ 



+ C F \n — In — - C F ln 2 ^ - ^C F ln ^ 



14 




2t 



2T, 



<-m 



-u 



S J \ s J \ s 
^ln ( tu \ (^ 2 + 2 ^) ln 2 M 




2s 2 




(m 2 + 2ms) 
2? 



u 



In 2 — 



In the above expressions for c 99 and c 99 , T qg and T qq are given by 



(3.11) 



-Q9 



qq 



= {s 2 + t 2 )/s 2 
= (t 2 + u 2 )/s 2 . 



(3.12) 



Note that, at NNLO, all the NNLL terms, i.e. O ([ln(s 4 /p 2 .)/s 4 ] + j terms, 
have been derived by matching with the exact NLO cross section. The scale- 
dependent terms at that accuracy can also be derived by the straightforward 
expansion of the resummed cross section. The result is the same, which 
provides a nice cross check. As an additional check, an explicit variation 
of the hadronic cross section with respect to scale yields zero up to (1/54)+ 
terms. From the expansion of the resummed cross section we can also derive 
the (1/54)+ terms containing squared logarithms of the scale. These terms 
in the coefficient c\{qg) are 



In' 



c 




t)+^(v) + f 



(3.13) 



while the corresponding terms in the coefficient c 2 (qq) are 




(3.14) 

Since we don't reproduce the full c 2 coefficients, we don't include them in 
our numerical results in the next section. 

Finally, we note that all gauge dependence cancels explicitly in the NNLO 
expansions, as expected. 
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4 Numerical results 

In order to examine the phenomenological consequences of the NNLO ex- 
pressions presented in this paper, the 0(aa 2 s ) program of Ref. |23| is used 
as a starting point. This includes the lowest order 2 — ^ 2 subprocesses and 
their virtual corrections, the 0(aa 2 s ) 2 — > 3 subprocesses, and a fragmenta- 
tion contribution calculated by convoluting 2 — > 2 parton-parton scattering 
subprocesses with appropriate photon fragmentation functions. The NNLO 
expressions are then added to the output of this program. It is important 
to note that the fragmentation contribution has not been modified by the 
inclusion of any resummed terms. On the other hand, for fixed target ener- 
gies the fragmentation component is generally small in the pt region spanned 
by the data. The results discussed in this section have been obtained with 
the renormalization and factorization scales set equal to each other. In each 
case two scale choices have been used, Pt/2 and 2px- The CTEQ5M [p4| 



parton distributions have been used. In each of the comparisons presented 
below, the calculation for the invariant cross section has been averaged over 
the appropriate region of rapidity or transverse momentum as specified for 
the data set being shown. 

The curves in Fig. 3 show the NLO and NNLO predictions for the in- 
variant cross section for direct photon production using a proton beam on a 
beryllium target. The predictions are compared with the experimental results 
at Pbeam = 530 GeV/c from the E-706 Collaboration P5]| . The theoretical 
curves have been multipled by a nuclear correction factor of 1.09 [p5[| . Com- 
paring the band formed by the solid curves with that formed by the dashed 
curves shows that the inclusion of the NNLO contributions reduces the scale 
dependence by about a factor of two. A similar reduction has been noted 
in [0. The origin of this reduction has been discussed in both [[F] and [p6 



Furthermore, over the px range shown in Fig. 3 the NNLO corrections for 
/i = Pt/2 are rather small. Thus, the inclusion of the threshold resummation 
corrections is unable to bring the predictions into agreement with these data. 
As noted in the introduction, it has been argued that initial state radiation 
can lead to non-zero values for the transverse momentum of the colliding 
partons and that the inclusion of such effects may bring the prediction into 
agreement with the data. However, to date this has only been treated in 
a phenomenological manner [|J and a more definitive theoretical treatment 
is needed before firm conclusions can be drawn. Some recent work in this 
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Figure 3. NLO and NNLO results for direct photon production in hadronic 
collisions compared to data from the E-706 Collaboration |25[ at Pbcam=530 
GeV/c. 



direction can be found in 27 



A similar pattern of behavior is shown in Fig. 4 for the 800 GeV/c 
results from the E-706 Collaboration. Again, the NNLO corrections are not 
sufficient to bring the predictions into agreement with the data. However, 
the reduction in the overall scale dependence remains. 

The UA-6 Collaboration has published data |28| for direct photon produc- 
tion in both pp and pp collisions. Results are available for both the rapidity 
and transverse momentum dependence of the invariant cross section. In Fig. 
5 the results for the pp collisions are shown. In this case the band formed 
by the NLO predictions overlaps the data at all values of px shown. The 
reduction in scale dependence provided by the inclusion of the NNLO terms 
is also evident. 

In Fig. 6 the corresponding results for pp collisions are compared with 
the UA-6 data. One can see that the predictions with /x = (with or 

without the NNLO terms) are compatible with the data and the now familiar 
reduction in scale dependence is evident. 

The resummation formalism presented in this paper allows one to study 
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Figure 4. NLO and NNLO results for direct photon production in hadronic 
collisions compared to data from the E-706 Collaboration PS| at Pboam=800 
GeV/c. 
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Figure 7. NLO and NNLO results for the rapidity distribution for direct 
photon production in hadronic collisions compared to pp data from the UA-6 
Collaboration E? 
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Figure 8. NLO and NNLO results for the rapidity distribution for direct 
photon production in hadronic collisions compared to pp data from the UA-6 
Collaboration fM. 



the dependence on both the transverse momentum and the rapidity of the 
photon. In Figs. 7 and 8 the predictions are compared to the rapidity dis- 
tributions published by the UA-6 Collaboration |28| for pp and pp collisions, 
respectively. Comparing the NLO and NNLO results shows that the rapidity 
dependence of the relative size of the NNLO corrections is mild. In each case 
the reduced scale dependence is evident. 



5 Conclusion 

In this paper the application of threshold resummation to direct photon 
production has been studied using the formalism of Refs. || [J. The expres- 
sions for the moments of the resummed cross section have been expanded to 
O (act's) and expressions for the leading, next-to- leading, and next-to-next- 
to-leading logarithmic contributions have been given. The NNLL estimates 
for the 0(aa s ) contributions have been incorporated into an existing next- 
to- leading logarithm 0(aa s ) program and compared to high statistics fixed 
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target data. The scale dependence of the results is markedly decreased by 
the inclusion of the NNLO terms. Furthermore, it was found that the size 
of the corrections is small when the factorization and renormalization scales 
are chosen to be /i = Px/2- Nevertheless, even with the inclusion of these 
threshold corrections, it is still not possible to achieve a good description of 
all of the fixed target data. Additional study of this process will be required, 
both theoretically and experimentally. 
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A Alternative form of the NNLO expansion 



In this Appendix we rewrite our NLO and NNLO expansions in terms of the 
variables v — (s + t)/s and w = —uj (s + t) used in Ref. 



The factorized form of the cross section for direct photon production can 
be rewritten as: 



E, 



d 3 a h 



d?p 1 



Yl / dx * dx b(Pfa/h A (Xa, Hi) 4>f b /h B (Xb, Hi) 
f 

w 1 da fa f b ^ lXi 



ttsv dv dw 



-(w,s,v,fi F ,a s (fi R )) . (A.l) 



The threshold region is given by w = 1. In general, a includes distributions 
with respect to 1 — w at nth order in a s of the type 



ln m (l 



w 



1 — w 



m < 2n — 1 . 



(A.2) 



Before presenting the NNLO expansions of the resummed cross section in 
terms of the variables v and w, we want to emphasize that crossing t and u 
mixes singular and non-singular terms in 1 — w . The non-singular terms are 
required in order to have the correct symmetrization under t u away from 
threshold. Therefore in producing numerical results we use the expansions 
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in terms of the variable S4 as presented in Section 3, where the symmetry is 
manifest and there is no mixing between singular and non-singular terms. 
It is convenient to express the results in a form similar to that in Eq. 
HI): 



VW{1 — V)S f—f = (TaAS, V) 



2\ \ k 



2k-l 



, fln l (l-wY 



1+2 





1 - w 



(A.3) 



The expressions for the two Born terms are 

1 



d fg( s ' v ) = ^ aa A T i9 v ( A ' 4 ) 



and 



Zqq(s, v) = —>Kaoi s e q T qq - (A.5) 



where T qg = 1 + (1 - v) 2 and T qq = v 2 + (1 - v) 2 . 
The coefficients for the NLO terms are as follows: 

c\{qg) = C F + 2C A 

c\{qg) = C F ^ + \nv^ - C A ln ^^^-(Cf + Ca) In 



(A.6) 

and 

~c\{qq) = (4C F -C A ) 

c\{qq) = -2C F In (^^) + C A ln(l - v) - ^ - 2C F In U^j 
c\m = C F (4 + m f^l) J^) + ^m. (A.7) 



2 V v 
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These results are in full agreement with the exact NLO results in Ref. [|1S 
Note that the expansion reproduces only the scale-dependent terms in c\. 
The full 5(1 — w) terms can be simply read off from the exact NLO calculation. 

The coefficients for the NNLO corrections for the qg — > ^yq subprocess 

are: 



cl(gg) 



c 3 (qg) 



^(C F + 2C A ) 2 
l(C F + 2C A ) 



C F -t + lnt> 



(C F + 2C A 
3 



- (C F + C A ) In 



cfln(^)+cf 



C A \n 



1 - V 

V 

Po(Cf 
2 V 4 



C F {-- + kic 
+ (C F + 2C A )^\n 



C A In 



1 -v 



C f 



(C F + C A )ln(^ 



+ -K(C F + 2C A ) 



+ A) 



c f 



4 



lnt> ) + 



In 



1-v 



71 



^-(C F + 2C A ) 2 
o 



(A.8) 



with 



and 



-19 



0i 



, l C F + C A ln 



1-v 



(A.9) 



-.19 



4 



7 
-( 

4 



+ 



Nr. 



T, 
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+ 



cv 

T, 



19 I 



--(v 2 - 2(1 - v)) \n 2 v + -(l- v) ln(l - v) - (1 - v) Inv 

+ -tt 2 (1 - v 2 ) + -(1 - v 2 ) ln 2 (l - u) - -(1 - v 2 ) lnv ln(l - v) 
4 4 2 

1 7T 2 

-(1 + 2«) ln(l -v) + —(1 - 4w + 5v 2 ) 

2 6 
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-(3v 2 - 2(1 - v)) lnv + - (3v 2 + 2(1 - v)) In 2 v 



+ - (v 2 + (1 - v) 2 ) ln 2 (l - v) - (v 2 + (1 - v) 2 ) lnv ln(l - v) . 

(A.10) 

The corresponding NNLO terms for the subprocess qq —>■ jg are: 



cl{qq) 
cl(qq) 



with I! 



and 



-(AC F -C A ) 2 
l(AC F -C A ) 



-2C F In 



1 - v 



2C F \n 



+ C A ln(l-v) 



+ 



{AC F - C A ) 



-2C F ln 



cf In 
1 - v 



+ cf 



00 



+ C A ln(l - v )-^--2C F hx[^- 



' 2 \ 1 2 



+ A) 



CVhi 



4 



(Inv - 21n(l -v)) 



7T 



^-(4C F -C A ) 2 





c a - Up 



3 , 

h In 

2 



1-v 



(A.ll) 
(A.12) 



0o 



/4 



it 



+ 



2T, 



H — - lnv — CV In v ln(l — v) 
6 

[v(2 + v) ln(l - v) + (1 - v)(3 - v) lnv 



+ (3v 2 + 2(1 - v)) In 2 v + (2v + 3(1 - v) 2 ) ln 2 (l - v) 
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2T, 



m \ (ll-16v(l-v)). 
d(1 - v) ln(l - v) + — In t> 



6 



+ (2, + 3(l-,» 



"(2 - i») 2 ' 
In u 



(A.13) 



Note that, at NNLO, we have derived all NNLL terms, i.e. all the 
O ([ln(l — w)/(l — w)] + ^j terms, by matching with the exact NLO cross sec- 
tion in Ref . ]13| . The scale-dependent terms at that accuracy have also been 
derived by the expansion of the resummed cross section, giving the same 
result. As an additional check, an explicit variation of the hadronic cross 
section with respect to scale yields zero up to [1/(1 — w)] + terms. From the 
expansion of the resummed cross section we can also derive the [1/(1 — w)] + 
terms containing squared logarithms of the scale. These terms in the coeffi- 
cient c\(qg) are 



In 2 



(C F + C A ) 



-C F -C A ln 



+ 



3A> 



(C F + C A ) In In (i 



(A. 14) 



while the corresponding terms in the coefficient c%(qq) are 



In 2 ( !%) C F 



C F 3-21n 



v J I 4 

^ (A. 15) 

Finally, we note that all gauge dependence cancels explicitly in the NNLO 
expansions, as expected. 
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